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Abstract 

An equation of motion for a one dimensional 
rocket is derived as a function of f(t), the mass 
flow rate into the acceleration chamber and v(x,t), 
the velocity distribution along the chamber, thereby 
including the transient flow changes in the chamber 
which have not been explicitly contained in other 
papers. The derivation of the mass density requires 
the introduction of the special time coordinate r, 
which satisfies the equation 0 = l-9r/9t-v9r/9x. The 
mass density expression is f(t-r)9T/9x. The 
equation of motion is derived from both classical 
force and momentum approaches and is shown to 
be consistent with the standard equation expressed 
in terms of flow parameters at the exit to the 
acceleration chamber. 

Introduction 

As propulsion systems become more advanced 
and sophisticated there is a tendency toward 
exercising more exact control over the mass flow 
rate into the engine and the acceleration imparted 
to the mass as it flows through the engine. Hence, 
the derivation of rocket equations of motion for 
advanced propulsion systems will have to be based 
on a rocket model that allows for the mass flow 
rate into the engine and the velocity distribution 
along the engine to vary with time. 

The objective of the analysis undertaken herein 
is the acceleration of the rocket expressed in the 
form : 

du/dt s F[M, x s ,x e (t), v(x, t). f(t>, x (t)] (1.1) 

where: 


in the literature is to treat the rocket as one entity, 
without separating the mass in the acceleration 
chamber from the rocket proper. 

The equations of motion which apply to a variable 
mass system have been a subject of debate among 
physicists, resulting in no less than eight 
articles!2-19 published in the American Journal of 
Physics during the period 1964-1966 alone. The 
author finds the following two approaches to be 
equivalent and applicable to the variable mass 
problem considered herein. 

1. The "momentum" approach: The time rate 

of change of momentum of all the particles 
involved must equal zero. 

2. The "force" approach: The sum of the forcqs 
acting on all the particles involved must equal 
zero. 

II. Derivation of the Mass Density 
in the Acceleration Chamber 

An analysis of the flow in the acceleration 
chamber (see Fig. 1) is required before the equation 
of motion can be obtained in the form of Eq. (1.1). 
The acceleration chamber in the rocket is used to 
accelerate the mass inserted into the chamber to 
the velocity at which it leaves the chamber. It is 
assumed that the mass enters the chamber with zero 
relative velocity. The velocity distribution along 
the chamber is 

dx/dt = v = v(x, t) (2.1) 

where vis a specified function of both distance along 
the chamber and time. The acceleration of the mass 
within the chamber can be obtained from this velocity 
distribution. 


M — Initial mass of the rocket. 

x s — Location of the inlet of the acceleration 

chamber. 

x e (t) — Location of the end of the acceleration 
chamber. 

v(x,t) — Velocity of the mass in the ac- 
celeration chamber with respect to the 
chamber. 

f(t) — Mass flow rate into the acceleration 
chamber. 

x C p(t) — Location of the center of mass of 
the rocket proper (the rocket minus the ac- 
celeration chamber). 


r 



4-fe 


AX 


Particle A' 






t = Time for particle A to move from x g to x 

At = Difference in times that particles V: b and d 
entered chamber 

Figure 1. Acceleration chamber notation 


The above quantities completely define the rocket 
herein considered and must be initially specified. 
The acceleration is derived for the time period 
before the acceleration chamber is full as well as 
for the period after mass has stopped entering the 
chamber. -The extensive literature 1 which 
considers the variable mass flow-variable exhaust 
velocity rocket does not contain an equation of motion 
in the form of Eq. (1.1), which would be the 
parameters normally specified. The approach taken 


The amount of time t that the particle at point 
x q at time tg has been in the acceleration chamber 
can be found by numerically integrating the equation 


dt/dx = 1 / v 


( 2 . 2 ) 


between the limits of x=xq, t=tg and x=x s , t=tg- 
r is then obtained from 


(2.3) 


1 



If the variables in Eq.(2.1) are separable, the 
following equation defines r 

l o x o 

/ .dt , / f- (2.4) 

V r 

where a and b are defined by 

v = ab = a(t) b(x) (2.5) 

As xq and tp are completely arbitrary, the general 
expression for r is 


As time elapses, this particle moves down the 
chamber. t increases uniformly with t. Hence, 
the total differential of T equals zero. 

dT = 0=-§f dx + ^dt (2.15) 

Substituting Eq. (2.14) into Eq. (2.15) gives 

o = 1 - ^7 - v (2.16) 

As Eq. (2.13) was for an arbitrary x, the partial 
differential of Eq. (2.16) with respect toxmust equal 
zero. 


r= r(x - “ ( 2 . 6 ) 

The mass flow rate into the acceleration chamber 
is given by f, a specified function of time. This 
function, along with the velocity distribution, 
determines the mass density in the chamber. 
Consider the mass Am that occupies a small length 
Ax centered at x. An expression for this mass is 
obtained by multiplying the mass flow at the time 
that the particle at x left the inlet of the chamber 
by an interval of time Ar. 


Am = f*Ar (2.7) 

The functions f and df/dt when starred should be 
modified by substituting t- t for t. For exanjple, if 
t=t, df/dt=3t, then f*=(t-r)' :i , (df/dt)*=3(t-r) ? . The 
time interval Ar is the difference in times that the 
particles on the right and left sides of Ax entered 
the inlet of the acceleration chamber and is given 
by (see Fig. 1) 

( 2 . 8 ) 


The mass density p is defined as the mass per unit 
length, yielding 


p- dm = lim Am . f * dr 
dx ~ Ax -♦ 0 Ax dx 


(2.9) 


Equation (2.9) may be demonstrated to be valid 
by using the equation of continuity 


?£+ o 

et ax 


Differentiating Eq. (2.9): 


Sp 

at 


dr 
dx at 


+ (1 - 


ar (§) 


at ' 8x dt 


dp £l_ ,ar.2 ( dL* 

85 g x 2-<3x» <dt> 


( 2 . 10 ) 


( 2 . 11 ) 


Inserting Eq. (2.11) into the equation of continuity 
gives 


0 = (1 “if“ v f5 ) 


dr / dfx* , i dr 
5x dt' 8x 8t 


3v dr 

' 55 "55 


;) f 


( 2 . 12 ) 


The easiest way to show that Eq.(2.12) is valid 
is to introduce a variable T. 


T 


- t - r 


(2.13) 


T is assumed to be associated with the particle that 
occupies x at time t. Partial differentiation of T 
results in 


ST = . Sr JjX , dr 
3x 9x ’ 8t _ 1 * 8t 


(2.14) 


_ ° T 
3x 9t 


9v dr 
dx 9x 



(2.17) 


Considering Eq. (2.16) and Eq. (2.17) it is apparent 
that Eq. (2.12) is valid. 

III. Derivation of Equation of Motion When 
Acceleration Chamber is Full 

This section contains a derivation of the equation 
of motion for the variable mass flow-variable 
exhaust velocity rocket for the time period when 
the chamber is full. Figure 2 shows the rocket 
nomenclature. All distances measured relative to 
the rocket are measured from a frontal position 
so that an x coordinate will always be positive in a 
rearward direction. 



Cp - Position of center of mass of the rocket 
proper 

c r - Position of center of mass of the whole 
rocket 

u,u cr “ Absolute velocities 
v, v e - Velocities relative to rocket 


Figure 2. Rocket nomenclature 


The mass of the rocket m r consists of the mass 
in the rocket proper plus the mass in the ac- 
celeration chamber 


m, = m p + j p dx (3.1) 

x s 

where x s and x e denote the start and end of the 
acceleration chamber. The distance x e is assumed 
to vary with time. Differentiating Eq, (3.1) using 
Leibniz's rule^O and considering x s to be constant 
results in 


dm r dm p dx e * 5, 

“dt dt + p e dt + J 8t 


(3.2) 


where p e is the density at x e . Using the equation 
of continuity (Eq. (2.10)) 


2 




dx 



(3.3: 


where v e is the velocity of the mass in the ac- 
celeration chamber at Xgand PgV s has been set equal 
to -dm p/dt. Combining Eq. (3/2) and Eq. (3.3) gives 
the ma H ss flow rate leaving the rocket 


Simplifying the above equation using Eqs. (3.1) and 
(3.3) gives 


2 

, d m x 

0= E-£E 

r dt dt 2 


J- 


/ %f 


dx 


(3.10) 


The mass of the rocket equals the mass of the rocket 
at time t g minus what has left the rocket 


dm 



dx ~ 


■ v > 


(3.4) 


The momentum p of the system (ie: the 

momentum of all the particles that were initially 
in the rocket) can be expressed as the sum of the 
momentum of the particles in the rocket proper plus 
the particles in the acceleration chamber plus the 
particles that have already left the rocket. The 
momentum shall be expressed in an inertial 
coordinate system which was at rest with respect 
to the rocket before it started to move. 


The rocket proper, denoted with subscript p, 
consists of a group of particles which are loosing 
mass at point x s . The center of mass x c of the 
rocket proper varies with time and is assumed to 
be specified. Appendix A contains a derivation 
(Eq.(A.9)) of the momentum of this body 


P CP , .. 

dt + X : 


dm 

P 

.ST- 


(3.5) 


where u is the velocity of a point fixed to the 
acceleration chamber of therocket. Using themass 
density expression (Eq.(2.9)), the momentum of the 
mass in the acceleration chamber is 


- J p (u - v) dx 


(3.6) 


Using the mass flow rate (Eq. (3.4)) out of therocket, 
the momentum p e of the mass which has left the 
rocket can be expressed as 

t 

P e = ) 'e (v e-^ )(U - V e )dt (3.7) 


t 

r dx 

m r = M - / p e (v e - -gjS.) dt (3.11) 


The mass of the rocket proper equals 


m = M - f f dt (3.12) 

P J 

0 

The following equation reduces to the continuity 
equation when expanded 


dfiv . I dv 
dt p y 3t 


+ 



&PV 2 

3x 


, dv 
dt 


3 pv 2 
3x 


(3.13) 


Combining Eqs. (3.10)-(3.13) and (2.9) yieldsthe 
equation of motion for the variable mass flow- 
variable exhaust velocity rocket 



f + v^)dx 
(3. 14) 

where fe and O r/9x) e are computed at x=x e . This 
is the nearest one can come to deriving a general 
expression in the form of Eq.(l.l). The only quantity 
in the above equation which is not contained in Eq. 
(1.1) is r, which can be explicitly obtained only for 
certain functions of v (see Section VIII for an 
example). 




t 


S 


where t s is the time at which mass first starts to 
leave the rocket. Combining Eqs. (3.5), (3.6) and 
(3.7) results in the total momentum of the original 
mass of the rocket. 


dm 


P = um 


-Sr 3 ^, 


dm 

-.sr 2 



t 


(u - v) dx + 


/ 


J dt 


(3.8) 


As the total momentum of the system is 
conserved, Eq. (3.8) may be differentiated and set 
equal to zero. The result is 


The above analysis is consistent with the results 
obtained by those who treat the variable mass 
flow- variable exhaust velocity rocket as a variable 
mass problem in which the center of mass of the 
whole rocket plays the leading role. The main result 
of these approaches is equation B-10, which can 
be expressed in terms of the nomenclature of Fig. 
2 . 

du d 2 m 

0 = m r ST + (x e " x cr> 2~ 

dt 

dx dx dm (3.15) 

+ < v + — ~ - 2 — — ) 

fit * A* ' A* 


du 

p dt 


dm 

dt 


dt 

’m dx 

if - + 'e <u - v e> ST 




+ (u - V ) 


where m r =m, u cr =u c , x cr =x c , x e =x 1 , and v e =vj. 
du cr /dt is the acceleration of the center of mass 
of the rocket in an inertial coordinate system and 
x cr is the position of the center of mass of therocket 
(see Fig.2). 


e c a e a 

du f , f , f <W . 

+ dT J v dx + » J Si dx - / dx 


It will now be shown that Eq. (3.10) is equivalent 
(3.9) to Eq.(3.15). The center of mass of the rocket is 
obtained from 
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m x = m x + / xodx 

r cr p cp J p 


Differentiating results in 


in the rocket proper and an integral of the mass in 
the acceleration chamber results in 

(3.16) 


du 

dt 


/ ' 4r dx - / 


dv 
P p dt 


f dv 

•’ J 'pit 


(3.24) 


p c p. 


dx f 

/ 


dp 

~5t 


in terms of the notation in Fig. 1 and where /> p is 
(3.17) the density in the rocket proper (for 0<x<x e ). 
Appendix A contains a derivation of the equation 


The following identity can be established by 
integration by parts using the continuity equation. 


. m. riv - 
■ at dx ■ ■ 


pv dx 


(3.18) 


Combining Eqs. (3.4), (3.17) and (3.18) and setting 
dm p /dt=- p g v g gives 


dm x dm x 
r cr _ p cp 


dm 

r 

“31 


dm 

-at® 


(3.19) 


dm 






dm 

. -at®- + / 



(3.25) 


Differentiating, setting P n „v„ dx /dt = 0, and using 
Eq. (3.13) gives Pe e 


x s + / p p tI dx + / 


p dt 


dx (3.26) 


Differentiating the above equation results in 


d 2 m x d 2 m x 

r cr . - 


* ^S- . x. 


~ nr -ar 


1/ 


dx 


(3.20) 


V ^ ♦ f 

e dt J 


dpv 

ST dx 


Combining Eqs. (3.4), (3.10) and (3.20) gives 


o . du A 

0 _ ~7u “ m — • k 

r dt r .. 2 


cr 


dx dx dm 
+ (v. +~ - 2—~- ) 


dt‘ 

cr 

dt ' dt 


<x = - x ) ■ 


(3.21) 


where it has been assumed that PpV^ is equal to 
zero atx equal to 0, x s and x ft . Equation (3.24) can 
be obtained by combining Eqs. (3.10), (3.13), and 
(3.26). 

IV. Derivation of Equation of Motion Before Mass 
Leaves the Acceleration Chamber 

The equation of motion for the time interval 
before the mass starts to leave from the acceleration 
chamber is found in a similar manner as was the 
main equation of motion (Eq.(3.14)). Themomentum 
expression can be obtained from Eq. (3.8) after 
eliminating the term representing the expelled mass 


The inertial acceleration (du cr /dt) of the center of 
mass of the rocket is 


du . 

Cr _ du 
dt " dt 


(3.22) 


Combining Eqs. (3.21) and (3.22) yields the desired 
result, Eq. (3.15). 

Equation (3.10) is also consistent with the ap- 
proach taken in references 5 and 6. Miele applies 
Newton's Second Law of Motion to the rocket 
regarded as a continuum in the following manner, 

*■* / (f - »,$- / f d ” 0.23) 

m ' ' m 


where F is the external force, du/dt is the ac- 
celeration of the origin (point fixed on the rocket), 
dv/dt is the relative acceleration with respect to 
the origin and the integral is to be summed over 
the mass of the rocket. Although Miele extends 
his approach further, it is easiest to derive Eq. 
(3.10) directly from Eq. (3.23). 


dm x 

p = um - P C P 

p eft 




(u - v) dx 


x 


s 


(4.1) 


where xj is the point reached by the first particle 
that entered the acceleration chamber. Xj is 
obtained by integrating Eq.(2.1). The mass of the 
rocket proper is given by Eq.(3.12). Differentiating 
Eq.(4.1), setting equal to zero and using Eq.(2.10) 
and the derivative of Eq.(3.12) results in 


0 = M 


x df 

x s Ht " 


w- s 


dpv 


(4.2) 


where vj is the velocity at xj. Combining Eqs.(3.12) 
(3.13) and (4.2) gives 


0 - M ^ ♦ (X. 


cp 


df 

V ht 


dx 

— 

dt 


- (M - J fdt ) 
0 


_C£ 


dr / 3v 9v . . / » q \ 

f 3x ( '5T dx ( 4 » 3 > 


Setting the external force equal to zero and 
separating the integral into an integral of the mass 
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The only difference between Eq. (3.14) and the 
above equation besides the obvious ones of 
interchanging xi for x g and the dropping of the 
integral in the first term is that in the above equation 
the derivative of xj with respect to time is vj, 
whereas in Eq. (3.14) x e is a specified function of 
time. 

V. Derivation of Equa t ion of Motion After 
Mass Stops Entering The Acceleration Chamber 

Once the mass stops entering the acceleration 
chamber, the mass remaining in the rocket proper 
remains constant. If t^ represents the time that 
the last particle entered the chamber, at time t it 
would have reached point X 2 found by integrating 
Eq. (2.1). Theintegral representing the momentum 
in the chamber need cover only the range from X 2 
to x e . Equation (3.8) for this case becomes 

X t 

dx f G r dX 

P =m p ( u--d T R) + J Plu-vldx+J P e <v e - -gj-) <U - v e > dt 
x 2 l s 

(5. 1) 

The mass of the rocket can be expressed 

X 

e 

m r =m p + / ' dx (5.2) 

x 2 

Differentiating Eq. (5.1), setting equal to zero and 
using Eqs. (2.10) and (5.2) results in 


where: m r = M if no mass has left the rocket. 
Otherwise 

t 

a dx o 

™ r = M -i f e <Jj>e (v e--dT >dt (6.2) 


x a equals x if mass is currently entering 
the acceleration chamber. Otherwise it is 
the point occupied by the mass in the chamber 
closest to the inlet. 

xb equals x e if mass is currently leaving the 
acceleration chamber. Otherwise it is the 
point occupied by the mass in the chamber 
closest to the exit. 

t a equals t if t< t^; otherwise t a = t^. 

The points x a and Xb when not equal to x s or x e 
respectively, are functions of time; they move with 
a velocity equal to the velocity in the acceleration 
chamber at that point. 

VII. Special Cases 

Several rocket models lead to simplified density 
expressions and corresponding simplified equations 
of motion. In the following examples both x cp and 
x e are assumed to be independent of time and x s 
is assume to occur at x cp . In addition, the reference 
point in the rocket is assumed to coincide with x s . 
Attention is restricted to the full acceleration 
chamber situation. 


. du ^ 

0 = m r -ar - m , 


d 2 x c „ 


x e A. Exhaust Velocity that Allows Separation of 

! dx Variables 

1 at (5.3) 

x 2 Assume that v is given by 


Combining Eqs. (3.11), (3.12), (3.13) and (5.3) gives 


v = ab = a(t) b(x) (7.1) 

Differentiating Eq. (2.4) and using Eq. (2.9) gives 




< v e--df> dt 


du 

dt 


- (M - j f dt) - / 


■ dr . 3v 9v,..j 

Si at + v ai ) dx 


(5.4) 


The differences between Eqs. (3.14) and (5.4) are 
the changes in the limits on the integrals and the 
dropping of the terms containing x s and f. 

VI. General Equation Of Motion 


3r/3x = 1/a b, fi = f /ab 


(7.2) 


Equation (3.14) reduces to 

X 

f 6 * 

du / f ,da , 2 db, . (7 

-at - J -« <dt + a s’ dx 

a 


Assuming further that f=a, the density P - 1/b and 
Eq. (7.3) reduces to 



0 = (M ■ 


, dt) ■ 


da 

dt 


(7.4) 


By comparing the equations of motion (Eqs. 
(3.14), (4.4) and (5.4)) for the different time periods 
considered in the previous three sections, the -fol- 
lowing general equation of motion for the variable 
mass flow-variable exhaust velocity rocket can be 
obtained 


0 = m 


du + (x 

■at <x « 


cp 


dx 

■ x s> £ + 2f 


d 2 X 


- (M - / f dt) 


T - I f si < sr + v 5i , dx 

dt 


( 6 . 1 ) 


B. Constant Exhaust Velocity Rocket 

The constant exhaust velocity rocket implies that 
the exhaust velocity doesnot change with timewhile 
the mass flow is variable. Hence, 


v = b = b(x) (7.5) 

Differentiating Eq. (2.4) with a equal to one gives 

3r/a x = i/b (7.6) 
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Eqs. (2.9) and (3.14) reduce to 


o ■ <M - / f* *) f- - / 

t o 


where r= f dx_ 

n b 


C. Constant Mass Flow Rate Rocket 

The constant mass flow rate rocket implies that 
the mass flow into the acceleration chamber does 
not change with time while the exhaust velocity does 
change with time. Hence, 


where c is a constant. Eqs. (2.9) and (3.14) reduce 
to 


t , x e 

M - c / <l5>e v e dt ] + 


D . Steady State Rocket 

Both the mass flow rate and exhaust velocity are 
independent of time for the steady state rocket. If 


f “ c . V = b = b(x) 

then Eqs. (2.9) and (3.14) reduce to 


(7.10) 


r = t - [t 2 - 2x 1_n / (l-o>] 1 

Eliminating x between Eqs. (8.2) and (8.3) gives 

, , r „ , n 2n-l 

+ nt 2 ij " [ t 2. (t _ r ,2] mT ( 8 . 4 ) 


Hence, to avoid infinite accelerations at the inlet 
(i.e.: when t goes to zero), n must lie in the range 
l/2<n<l. 

Differentiating Eq. (8.3) yields 

f = x - n ( t 2 . 1 i _ x 1 - n ,- 1/ 2 (8 - 5) 

The density (Eq. (2.9) and mass flow rate are 

,-./n (8.6) 

p ~ 1/X ; fiv s t 

Equation (4.4) reduces to 

° = M |-vV> 2 (8,7) 

Integrating the velocity to get the time-position 

history of the first particle that enters the chamber 

results in 

Xj « £(l-n) t 2 /2j 1/2 (8.8) 

Inserting into Eq. (8.7) and integrating gives 





3-n r 

n ._ J_ 1 

M - c (t - t ) 
s 

du , 

dt ‘ cb e 

(7.11) 

„ _ i 1-n t 1-n 
M 3-n 1 1 

1 - n + (ifa) '- n J 


where b Q is based on x = x . Integration of 
Tsiolkovskii's equation gives e 

u - u s = b e lo g {M/(M-c(t-t s )]} (7.12) 

where u g is the velocity at time t g . 

VIII. Numerical Example 

This section contains an example of when the 
variables in v are separable and such that Eq. (2.4) 
can be integrated. The same assumptions 
concerning x C p and x are made as was assumed 
in Section VII\ Let the mass flow rate into the 
acceleration chamber and the velocity be specified 
by 

f = t ’ V = txn (8.1) 

The acceleration of the mass in the acceleration 
chamber is given by 

^ = ir + ^ = x" + nt 2 x 2n - 1 (8.2) 

Carrying out the integration of Eq. (2.4) resultsinr 


This equation is valid until the chamber is full at 
time t_. This time is found by solving Eq. (8.8) 
for t after substituting x e for xj. 


t s * [2x e I_n / (l-n)] 1/2 (8- 

For times in excess of t g , Eq. (3.14) reduces to 


( 8 . 10 ) 


( 8 . 11 ) 


0 = (M - t 2 /2 + t 2 /2) - x - x n t 2 

s at e e 


Integrating 

“ - u s ■ Xg (a - b)/c + 2x e n (-t + ac/2 + t s - bc/2) 

( 8 . 12 ) 

where 

c = (2M + tg 2 ) 1/2 , a = log [(c + t)/(c - t) ] 
b = log [(c + t g )/(c - t g )] 

This equation is valid until mass stops entering 
the chamber at time tj. For greater times the 
equation of motion is found using Eq.(5.4) 
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0 - <M - t 2 /2 + t 2 /2) % - x" t 2 - x + x" t 2 + x„ 


(8.13) 


where X£ is obtained by integrating the velocity. 


[(l-n) (t 2 - 1 , 2 ) n ] 


l/(l-n) 


(8.14) 


3. Snyder, V., "Comment on "Motion of the Center 
of Gravity of a Variable-Mass Body", "AIAA J.3, 
571-572 (1965). 

4. Thomson, W., "Equations of Motion for the 
Variable Mass System", AIAA J.4, 766-768 (1966). 

5. Miele, A., Flight Mechanics (Addison Wesley 
Published Co., Inc., Reading, Mass., 1962) pp 23-31. 


As Eq. (8.13) is not readily integrable, it can 6. Halfman, R., Dynamics, Vol. 1 (Addison Wesley 

be solved numerically on a computer. Figure 3 is Publishing Co., Inc., Reading, Mass., 1962) Section 

a plot of u verses t for the following conditions 4-3. 


7. Meriam, J., "Variable-Mass Dynamics", J. End. 
(8.15) Ed., 51, 241-243 (1960). 


m p at l s = • 75 ' m p at 4 = . 25 



TIME 

Figure 3. Velocity and distance versus time for 
rocket defined in Eq. (8. 15) 


The length of the acceleration chamber equals 
.00316 and is found by setting and t equal to x e 
and t g in Eq. (8.8). The time that the last particle 
leaves the acceleration chamber equals (/2and is 
determined by setting x 2 and t equal to x e and tg 
in Eq. (8.14). The distance the rocket has moved 
is found by numerically integrating the equation 
dy/dt = u and is also shown in Fig. 3. 
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Appendix A 

Deprivation of the Momentum of 
A Variable Mass Body In Terms 
Of Its Center Of Mass 


The derivation of the momentum of a variable 
mass body in terms of its center of mass can be 
based on a discontinuous medium ^ or on a continuous 
medium 1*9. To be consistant with the general 
approach taken in this paper, a body with a continuous 
medium will be used. The body is assumed to lose 
mass at both point x^ and x,, with x^ and x 2 moving 
with time (see Fig. 4). The mass leaves with the 
local velocity v of the mass in the medium. 


The location x c of the center of mass is given by 


m x 


c 



x p dx + 



x p dx 


(A.l) 
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Q 

>c 

^ V, 

— ►dx/dt 

_v 2 
f— ►d«2 

— 

J 

X 

dt 

1 -I 

1 1 


c - Position of center of mass 
u - Absolute velocity 
v - Velocity relative to body 


Figure 4. Variable mass body notation 

There is a discontinuity in the mass density at xi 
due to the mass leaving at that point. If /oqi and 
^12 refer to the mass densities on the left and right 
sides of xi, differentiating Eq. (A.l) results in 


x i *01 

dXj 

dt 

X 2 ' 

dx 2 
2 dt 


dx. 

x ! 

x 

x 3 1 

X 9 



L 

r 2 

dx + f 

■'v 

x If dx 


(A. 2) 


Using Eq. (3.18) the above equation can be reduced 
to 


dmx 
c 

dt 



+ x i < 'oi 


) 


f 0 ' vdx + L 


p v dx 


The mass in the body equals 


(A. 3) 


m 



P dx + 



dx 


(A. 4) 


Differentiating and using the continuity equation 
gives 


dm 

dt 




o) ( 


dXj 

“3T 


' v ,> 


(A. 5) 


The first term on the right hand side of the equation 
represents the mass flow leaving point 2 and the 
second term the mass flow leaving point 1 
Combining Eqs. (A.3) and (A.5) results in 


dx c dx ? 

m dt- + *2 <^r-v 2 ) <* c -x 2 ) 

dXj X 1 x 2 (A # 6) 

+ ( ~dT~' V l > ( 'oi " ' > 12 ) <x 2 ‘ x i* 3 f g e v dx + f p v dx 

X 1 

The momentum of the body is defined to be 



P - v) dx + 



(A. 7) 


where u is the velocity of a point fixed on the body 
Comparing Eqs. (A.6) and (A.7) gives 


l=u / pdx+u / 
0 \ 


p dx - m 


dt 


dx 2 1 dx 

+ "2 ( -ar - v 2> < X J - x c > + <-3T - Vj ) </> 0] - ^ 12 )(Xj - x c ) 


(A. 8) 


Equations (A.6) and (A. 8) can be applied to the 
rocket proper section of Fig. 1. Using the notation 


of Fig. 1 and Eq. (A.5) with Pq (dx 2 /dt-V 9 ) set equal 
to zero, Eqs. (A.6) and (A.8) reduce to 


dx s A e 

1 ~ar + <x c * x s ’ if = / * v dx + f p v < 


(A. 9) 


P = um - m - 


• + (x - xj 


dm 


dt ~s " A c' *dT 

Appendix B. 

Equation of Motion For A Variable Mass B ody 
In Terms Of Its Center Of Mass 


1 nis appendix contains a derivation of an equation 
which has been frequently mentioned in the 
literature as the equation of motion of a variable 
mass body (of which the rocket is one example). 
The derivation follows the approach taken in the 
main section of this paper (differentiating the 
momentum of the system and setting the resulting 
expression equal to zero). The body considered is 



— r * — — 

C r*» 

fl 

i » - dm 

U-, 1 

— ►dx, 

ST 

— x c — * 


1 *1 H 


c - Position of center of mass 
u - Absolute velocity (*dr/dt) 
v - Velocity relative to body 
dm/dt - Mass flow rate at x^ 


Figure 5. Variable mass body notation 
The mass within the body is 


Assuming that the body loses mass at point x-i which 
varies with time, Eq. (B.l) can be differentiated 
using Eq. (2.10), resulting in 


dm _ 


,- at— V 


(B.2) 


where v j is the velocity of themassbeing exhausted. 


The momentum of the body plus the momentum 
of the mass that has left the body equals 



dx i 

<v l "dT 1 (u - v j> dt 


(B.3) 


where u is the velocity of the body and t s is the 
time that the body started to lose mass 
Differentiating Eq. (B.3) using Eqs. (2.10) and (B.2) 
gives 


fr° = 


du 
1 dt 



dpv 

“3r 


dx 


(B.4) 


The location x c of the center of mass is given by 

x i 

mx c ~ f x p dx (B.5) 

^ J o 

Differentiating gives 


drrix c dx i P 1 a, 

' x i p i -ar + / Q x Jt dx 


(B.6) 


where dx c /dt is the velocity of the center of mass 
with respect to the body. Using Eqs. (3.18), (B.2) 
and (B.6) results in 
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dmx 


dm J 


(B.7) 


p v dx 


dt o 

Differentiating once more gives 

dx. 


. , dm , 1 / 

+ <-dt' + "i v i ) ~dT + J 


dpv 

eft 


(B.8) 


dt dt 0 

The acceleration of the center of mass point is 


c _ du 

dt " dt 


d x 
c 

m 2 


(B.9) 


Combining Eqs. (B.2), (B.4), (B.8) and (B.9) results 
in 

du .2 dx, dx 


0 = m 


dt 


, , , d^m , , , UA 1 0 - A c » dm 

+ ( - xj — ^ + ( v, + - 2 - 


dt 


dt 


~3r ' ~ar 


(B.10) 


Introducing the notation r c , r g and v (absolute 
velocity) with the following equations (see Fig. 5) 


i 

dx. 


dr 

dt 


dr dx 
e c 

dt * dt 


dr 

dt 


(B.ll) 


results in the equation of motion as usually specified, 

du„ dr dr 

+ (2 Tf-Tf- v >^ (B.12) 
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